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Abstract.  Variation  of  the  compre-sional  velocity  with  depth  in  the  C-region 
of  the  upper  mantle  in  western  North  America  is  studied.  Apparent  velocities 


of  first  arrivals  across  the  Tonto  Forest  Array  in  Arizona  jiatf^-btten  determined 
for  about  70  shallow  focus  earthquakes.  The  epicenters  range  from  10  to  30  degrees 
in  distance  and  are  mostly  south  of  the  array. 

The  method  gives  the  slope  of  the  travel  time  curve  directly,  the 
parameter  required  for  &  velocity-depth  determination.  For  this  distance  range 
it  is  found  that  although  the  slope  of  the  P  travel  time  curve  decreases  sub¬ 
stantially  with  distance,  i.e.,  increasing  apparent  velocity,  the  variation  is 
not  uniform.  Two  relatively  abrupt  changes  are  observed  at  distances  of  about 
17  and  24  degrees.  These  are  most  readily  interpreted  as  two  second  order  dis¬ 
continuities  in  the  region  of  the  mantle  above  1000  km. 

Assuming  an  average  crustal  structure  as  derived  from  seismic  refrac¬ 
tion  measurements  and  taking  the  uppermost  200  km  of  the  mantle  given  by  Gutenberg 
or  Jeffreys,  various  possibilities  are  explored  for  the  form  of  variation  of  the 
P  wave  velocity  with  depth.  The  discontinuities  in  the  rate  of  change  of  velocity 
with  depth  which  cause  the  observed  breaks  in  the  apparent  velocity  curve  are 
found  to  be  at  depths  of  about  320  and  640  km.  Although  no  absolute  times  are 
required  in  the  method  the  travel  times  for  the  various  derived  models  are  computed 


and  compared  with  standard  tables, 
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INTRODUCTION 

The  most  abundant  and  direct  source  of  data  regarding  the  structure 
of  the  earth’s  interior  is  the  travel  time  of  various  seismic  body  phases.  The 
Herglotz-Wiechert  method  used  in  the  interpretation  of  body  waves  requires  a 
knowledge  of  d t  /d  &  as  a  function  of  distance.  This  is  ordinarily  obtained 

from  the  slope  of  the  observed  travel  tim*  curve  and  involves  the  smoothing  of 

*• 

observations  made  at  seismic  observatories  with  different  crustal  structures. 

Application  of  this  method  has  yielded  a  number  of  velocity-depth 
profiles  which  are  similar  except  in  the  upper  1000  km  of  the  mantle  and  a  small 
region  of  the  transition  zone  of  the  core.  There  are  various  reasons  for  the 
differences  in  the  derived  structures  of  the  upper  mantle  including  real 
variations,  errors  of  apparent  velocity  introduced  by  smoothing  of  scattered 
data  and  disagreement  in  the  interpretation  of  data.  The  existence  or  absence 
of  low-velocity  layers  and  discontinuities  are  the  main  differences  between  the 
standard  models  of  Gutenberg,  Lehman,  and  Jeffreys.  These  are  also  the  complica¬ 
tions  which  lead  to  theoretical  difficulties  in  the  body  wave  travel  time  method. 

The  models  of  Gutenberg  and  Lehman  have  an  upper  mantle  low-velocity  region; 

Jeffrey’s  model  does  not.  The  models  of  Lehman  and  Jeffreys  have  a  strong  discon¬ 
tinuity  above  a  depth  of  400  km;  Gutenberg’s  model  is  smooth  throughout  the 
mantle.  The  scatter  of  travel  time  data  and  the  complicated  amplitude  pattern 
in  the  epicentral  region  of  15°  to  about  28°  is  partly  responsible  for  the  ambiguities 
of  interpretation. 

The  first  evidence  for  the  sudden  change  of  the  slope  of  the  travel 
time  curve  of  P  waves  near  the  epicentral  distance  of  20  degrees  was  presented 
by  Byerly  (  1926^.  In  a  later  analysis  of  the  Nevada  earthquake  of  December  20,  1932 
with  considerably  more  date  (Byerly,  1935)  he  noted  that  the  P  travel  times 
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showed  two  discontinuities  at  epicentral  distances  of  about  17  and  28  degrees, 
the  latter  being  clearer.  In  the  last  few  years,  measurements  have  revealed 
additional  evidence  for  discontinuities  in  the  upper  1000  km  of  the  mantle.  A 
brief  account  of  these  observations  is  summarized  by  Anderson  (1965^  . 

In  recent  years  a  number  of  seismometer  arrays  have  been  installed 
around  the  world,  some  having  linear  dimensions  of  several  kilometers.  The 
primary  purpose  of  these  installations  is  the  improvement  of  signal  to  noise 
ratio  to  increase  the  capability  of  the  detection  of  small  events.  Each  array 
is  pra  Lcally  equivalent  to  several  stations  with  essentially  identical  under¬ 
lying  crust  and  having  the  same  azimuth  from  the  epicenter.  They  can,  therefore, 
be  used  for  direct  measurement  of  d tT /d  &  of  various  seismic  phases. 

The  linear  extent  of  an  array  which  can  be  used  for  this  purpose  is 
restricted  on  one  hand  by  the  required  resolution  of  the  time  of  arrival  of  the 
signal,  and  on  the  other  hand  by  the  desired  homogeneity  of  the  underlying  crust 
and  upper  mantle.  While  the  latter  restriction  depends  on  the  degree  of  lateral 
variation  in  the  uppermost  50  to  100  km  of  the  earth’s  structure,  the  former 

varies  with  the  phase  under  consideration  and  with  epicentral  distance. 

6V 

In  this  paper  the  seismograms  of^shallow  earthquakes,  selected  from  among 
more  than  200,  and  one  underground  explosion,  recorded  at  Tonto  Forest  Seismological 
Observatory  (TFSO) ,  Arizona,  are  studied.  The  d^/dA  of  the  first  arrivals 
in  the  P  group  of  these  earthquakes  across  the  array  were  measured  and  examined 
for  possible  indications  of  discontinuities  in  the  region  of  the  mantle  between 
225  and  800  kilometers  depth. 

Details  of  the  installation  and  operation  at  TSFO  are  published  in 
the  semi-annual  progress  reports  of  the  observatory.  Data  are  read  from  the 
recordings  of  the  short  period  vertical  seismometers  installed  at  the  end  points 
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of  a  crossed  linear  array  (Z63,  Z67,  Z70  and  Z74  in  Figure  1).  The  array 
consists  of  two  perpendicular  lines,  each  of  which  is  approximately  10  km  lone 
and  contains  eleven  seismometers.  The  instruments  are  of  the  Jchnson-Matheson 
type,  running  with  unusually  high  magnification  of  650  K  at  1  cps. 

In  Figure  2,  the  recordings  of  first  P  arrivals  of  two  of  the  earth¬ 
quakes  are  reproduced.  The  first  earthquake  has  a  magnitude  of  5.4  and  is 
located  13.3  degrees  toward  the  northwest  of  the  array.  The  second  one  has  an 
e picentral  distance  of  20.3  degrees  to  the  southeast  of  TSFO  and  a  magnitude 
of  4.1.  The  clear  later  arrivals  are  an  additional  source  of  information  regarding 
structure  of  the  upper  mantle  which  we  do  not  use  in  the  present  paper  because 
of  uncertainties  in  the  depths  of  the  sources. 

ANALYSIS  OF  DATA 

A  list  of  the  data  is  given  in  Table  1.  All  but  the  numbers  listed 
in  the  last  column  are  reproduced  from  the  Geotech  monthly  bulletins,  and  are 
based  on  the  USCGS  preliminary  reports.  The  measured  slopes  are  given  in  the 
last  column.  Only  two  of  the  earthquakes  have  listed  focal  depthsof  over  50  km 
for  which  appropriate  corrections  in  the  epicentral  distances  are  made. 

The  distribution  of  the  epicenters  arcund  TSFO  is  exhibited  in  Figure  1. 
The  epicentral  distances  range  from  12  to  28  degrees.  About  twenty  per  cent 
of  the  azimu_hs  are  to  the  northwest;  the  remaining  80  per  cent  stretch  fron;  due 
south  to  the  southeast  of  TSFO. 

A  plot  of  the  measured  slopes  is  shown  in  Figure  3.  Solid  circles 
represent  earthquakes  located  to  the  northwest  of  TFSO.  The  azimuths  of  the  earth¬ 
quakes  shown  by  the  hollow  circles  range  from  due  south  to  the  southeast.  The 
point  marked  by  6  belongs  to  the  Salmon  underground  explosion  in  Mississippi 
(<~4=  18. 5°,  Az  =  93.7°).  For  this  explosion  the  slope  was  measured  on  the 


horizontal  components. 
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In  the  same  figure  the  tabulated  Jeffreys  slopes  for  shallow 

earthquakes  (h=0,00  R)  ,  computed  from  the  Jef freys-Bullen  tables  (l?58)  are 
smoothed  and  shewn  by  the  dashed  line.  It  is  noted  that  a  variation  of  slopes 
such  as  the  one  given  by  the  Jeffreys  model  is  not  in  close  agreement  with  our 
measured  values.  Our  data  fails  to  indicate  any  sudden  decrease  in  cl  T/A-k 
at  il  =  20°.  Instead,  there  are  indications  that  the  assumed  discontinuity 
at  this  distance  is  the  compounded  effect  of  the  gradual  decrease  of  the  slope 
in  the  spicentral  range  between  15  and  25  degrees  which  may  include  two  less 
pronounced  discontinuities  at  epicentral  distances  of  about  17  and  24  degrees. 

Evidence  for  the  existence  of  twin  discontinuities  with  close  corres¬ 
pondence  to  the  ones  suggested  here  have  previously  been  presented  by  a  number 
of  investigator:,  (see,  for  example^  Doyle  and  Webb,  1963).  With  the  commonly 
used  technique  of  few  observations  several  degrees  apart,  it  is  extremely  difficult 
to  see  the  fine  structure  of  the  travel  time  curve.  Furthermore,  when  the  ob¬ 
served  points  are  substituted  by  straight  lines  the  information  concerning  small 
abnormalities  are  lost. 

The  data  presented  in  Figure  3,  though  illustrating  the  general  trend, 

is  quite  scattered.  This  is  in  particular  true  for  distances  less  than  15 

degrees.  Most  of  the  earthquakes  in  this  range  are  from  off  the  coast  of 

northern  California  and  Oregon  (solid  circles).  Paths  of  these  earthquakes 

traverse  the  anomalous  region  underlying  the  Sierra  Mountains.  Also  some  scattering 

arises  from  the  fact  that  may  still  appear  as  the  first  arrival.  Miss  Lehman 

(1964)  observed  that  to  the  northeast  of  the  Gnome  explosion  was  a  first 

arrival  to  epicentral  distances  as  far  as  19  degrees.  The  P  velocity  also 

n 

varies  rather  rapidly  in  different  regions  of  the  western  United  States.  Large 
amplitude  later  arrivals  may  also  have  been  picked  as  first  arrivals. 
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Other  factors  contributing  to  the  scattering  are  errors  in  measure¬ 
ments  and  uncertainty  in  the  epicentral  distances.  By  analyzing  only  superior 
seismograms  with  low  background  noise  and  prominent  first  arrivals  we  attempted 
to  minimize  reading  errors.  Nevertheless,  such  errors  may  result  in  uncertainties 
as  large  as  1.0  sec/deg  in  the  measured  slopes.  More  accurate  location,  on 
the  other  hand,  requires  an  even  distribution  of  recording  st  tions  around  the 
epicenter.  This  is  not  realized  for  the  earthquakes  under  consideration,  because 
of  their  geographic  location.  A  study  of  this  nature,  therefore,  should  include 
a  large  number  of  observations  so  that  a  statistical  treatment  of  the  data  becomes 
possible. 

Two  alternate  methods  were  adopted  to  smooth  the  data  prior  to  interpre¬ 
tation.  In  the  first  method  all  the  observations  were  fitted  by  a  single  least 
square  polynomial  of  degree  9,  giving  a  standard  error  of  estimate  of  0.43  sec/deg. 

In  the  second  method  the  data  were  treated  in  three  groups  with. breaks  at 
distances  of  17  and  24  degrees.  Each  group  was  fitted  by  a  polynomial  of  degree  5. 

The  standard  error  of  estimate  for  the  whole  range  in  the  latter  case  is  reduced 
to  0.35  sec/deg.  The  discontinuous  solid  curve  in  Figure  3  resulted  from  the 
second  procedure. 

Due  to  the  uncertainty  in  the  nature  of  the  first  arrivals  and  high 
scattering  of  the  data  at  distances  less  than  15  degrees  we  made  use  of  two  well 
established  models  for  depths  corresponding  to  this  epicentral  range.  The  uppermost 
200  km  of  the  mantle  is  alternately  assumed  to  have  P  velocities  suggested  by 
Jeffreys  (see  Bullen, (1963)  and  Gutenberg  (1959)).  Based  on  the  preliminary 
results  of  recent  seismic  refraction  work  in  the  recording  region  (Warren  et  al,  1965), 
a  single  layered  crust  of  30  km  thickness,  having  an  average  compressional  velocity 
of  6.1  km/sec  is  also  assumed. 
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Under  these  assumptions  and  a  choice  of  continuous  or  discontinuous  treat¬ 
ment  of  the  data,  a  number  of  possible  velocity-depth  profiles  are  derived  for 
depths  up  to  800  km.  The  method  of  successive  stripping  is  applied  so  that  prior 
to  the  application  of  the  Herglotz-Wichert  formulas 

(r.  / r, 

the  parameters  involved  i.e.,  /A  ,  and  r(  are  properly  reduced  to  those  of 
the  stripped  earth.  Here  rQ  is  the  radius  of  the  earth,  p  parameter  of  the 
ray  (equivalent  to  the  slope  of  travel  time  curve),  v^  velocity  at  the  bottom 
of  the  ray  (r=r^)  emerging  at  an  epicentral  distance  A  ^  and  having  parameter 
Pj.  V  is  the  apparent  surface  velocity  at  distance  A  . 

All  the  computations  were  performed  cn  an  IBM  7094  computer,  and  results 
are  summarized  in  Table  2.  A  comparison  of  the  variation  of  P  velocity  in  the 
uppermost  800  km  of  the  mantle  derived  here  and  those  suggested  by  Jeffreys  and 
Gutenberg  is  illustrated  in  Figure  4.  Models  X  and  Y  are  constructed  on 
the  basis  of  the  hypothesis  of  gradual  decrease  of  d  t  /d  LS.  in  the.  epicentral 
range  of  15  to  30  degrees  and  differ  in  the  assumed  structures  between  M-dis- 
continuity  and  225  km  depth.  In  this  interval  X  coincides  with  the  Jeffreys 
model  while  Y  with  that  of  Gutenberg.  Model  Z  also  has  a  Gutenberg  upper 
mantle,  but  the  velocities  between  depths  of  225  and  800  km  are  based  on  a  d^/d  A 
curve  which  has  two  discontinuities  at  distances  of  17  and  24  degrees.  The 
depths  corresponding  to  these  discontinuities  are  close  to  320  and  640  km. 

DISCUSSION  OF  RESULTS 

As  is  apparent  in  Figure  4,  all  three  models  begin  as  the  continuation 
of  either  Jeffreys  or  Gutenberg  upper  mantle  (according  to  the  assumptions  made 
in  the  previous  section)  from  a  depth  of  about  225  km.  For  a  good  portion  of  the 
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profile,  models  X  and  Y  represent  an  average  to  th*  Jeffreys  and  Gutenberg 
velocities  and  together  with  model  Z  converge  to  a  value  which  is  higher  by 
about  0.2  km/sec  at  a  depth  of  800  km. 

In  the  sectional  smoothing  of  the  data,  one  point  at  the  epicentral 
distance  of  15.9  degrees^  for  which  the  observed  dt’  /d  A  is  12.9  sec/deg,  has 
been  assumed  to  be  a  second  arrival  and  is  thus  combined'  with  the  observations 
between  17  and  24  degrees  (see  Figure  3).  There  is  some  arbitrariness  in  the 
way  different  portions  of  the  curve  meet.  The  accuracy  of  the  present  data  is 
insufficient  for  deciding  whether  each  portion  continues  along  the  dash-dot  tail 
to  jo;  ...  the  adjacent  portion  smoothly,  or  there  is  some  overlapping.  Clear  later 
arrivals  on  some  of  the  records  suggest  that  the  travel-time  curve  is  multi-valued. 

For  the  sake  of  computational  convenience,  it  is  here  assumed  that  the 
neighboring  sections  of  the  discontinuous  curve  slightly  overlap.  The  corres¬ 
ponding  jumps  in  the  velocity  profile  are  quite  small  (0.03  and  0.06  km/sec  at 
320  and  640  km  depths,  respectively),  and,  in  the  absence  of  other  evidence  such 
as  reflections  from  such  first  order  discontinuities,  may  be  ignored.  In  sketching 
Figure  4,  we  have  taken  the  resulting  two  discontinuities  of  model  Z  to  be  of 
second  order. 

Any  suggested  earth  model  should  obviously  satisfy  the  observed  travel 
times.  The  three  models  resulting  from  this  study  are  based  on  measurements  of 
dt  /d&  ,  without  direct  reference  to  the  absolute  values  of  the  transmission 
time.  Indirectly,  however,  the  computed  distances  have  been  influenced  by  the 
selection  of  a  set  of  travel  time  curves.  The  bias  introduced  by  such  procedure 
is  not  expected  to  be  more  serious  than  other  sources  of  error  which  exist  in  the 
evaluation  of  distances.  Furthermore,  with  less  scattered  data,  it  is  possible 
to  remove  such  bias  by  a  process  of  successive  approximation. 

To  examine  the  compatibility  of  the  models  with  the  travel  times  which 
are  in  common  use,  we  have  computed  the  variation  of  the  transmission  time  with 
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distance  of  these  models  for  a  focal  depth  of  33  km  (see  Anderson  and  Julian. 

1965).  A  comparison  of  the  results  of  the  computed  time  of  first  arrivals  and 
JB  tables  is  made  in  Figure  5.  The  initial  values  of  the  residuals  (atZ.  =15°) 
can  be  explained  by  ,_he  difference  in  the  crustal  structures  and  also  that  between 
Gutenberg  and  Jeffreys  upper  mantles  (for  models  Y  and  Z) . 

It  can  be  seen  from  Figure  5  that  models  X  and  Z  are  in  better  agree¬ 
ment  with  JB  tables  than  model  Y.  The  residuals  for  model  Y  remain  positive 
throughout  the  distance  range  under  consideration  and  attain  a  na/imum  value  of 
3.3  seconds  atA=  23.5  degrees.  The  residuals  of  the  other  two  models  include 
both  positive  and  negative  values.  The  most  conspicuous  positive  range,  common 
in  all  three  models,  is  at  distances  between  about  21  to  28  degre_s.  Such  positive 
anomalies  based  on  various  assumptions  can  be  interpreted  to  indicate  that  for 
this  region  the  JB  travel  times  should  be  slightly  (at  least  one  and  possibly 
up  to  two  seconds)  increased  in  this  distance  range.  A  comparison  between  the 
residuals  of  models  Y  and  Z,  both  of  which  have  a  Gutenberg  upper  mantle, 
indicates  a  preference  for  a  discontinuous  d  1  /dA  variation,  such  as  the  one 

used  to  construct  model  Z,  over  a  continuous  one. 

M  o 

From  — =  28  on  the  residuals  of  X  and  Z  become  negative.  Such 
a  decrease,  in  all  three  models,  is  due  to  high  P  velocities  at  depths  near 
8o0  km  found  here. 


UNIQUENESS 

Although  we  feel  that  the  direct  measurement  and  analysis  of  apparent 
velocity  data  is  superior  to  the  usual  method  of  smoothing  and  differentiating 
travel-time  data,  particularly  for  the  elucidation  of  fine  structural  details, 
there  still  remains  a  great  deal  of  arbitrariness  in  the  interpretation.  In 
addition  to  the  complications  mentioned  previously  there  is  the  possibility  that 
the  three  regions  of  the  P  wave  arrival  data  which  we  infer  between  10  and  30 
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degrees  are  three  separate  branches  which  continue  into  the  adjacent  regions  as 
later  arrivals,  perhaps  even  coalescing  to  form  cusps.  Without  using  later 
arrivals  and  amplitude  data  it  is  impossible  to  discuss  the  sharpness  of  the 
discontii.ui-L^es .  While  the  gradient  of  velocity  in  each  depth  region  is  probably 
determined  correctly,  limited  only  by  the  scatter  of  the  data,  the  details  of  th* 
transition  regions  are  not  resolved. 

Our  method  of  analysis  can  be  considered  to  give  the  smoothest  variation 
of  velocity  consistent  with  our  tri-partite  d^  /d  A  curve.  The  model  CITll 
(Anderson  and  ToksBz,  1963  ;  Anderson.  1965)  also  has  two  major  discontinuities 
at  depths  close  to  those  determined  here.  Since  this  model  was  determined  from 
surface  waves  it  is  stepped  rather  than  smooth  and  represents  probably  the 
coarsest  variation  of  velocity  consistent  with  surface  wave  data.  The  actual 
nature  of  the  discontinuities  probably  lies  between  these  two  extremes.  The 
CITll  model  predicts  multi-valued  travel-time  curves  over  most  of  the  interval 
between  14  and  38  degrees.  Model  Z  has  only  minor  later  arriving  branches. 

Since  we  cannot  distinguish  between  these  two  kinds  of  models  on  the  basis  of 
first  arrivals  alone  we  content  ourselves  here  with  a  few  numerical  experiments, 
maiuly  to  study  the  effect  of  the  sharpness  of  the  discontinuity,  in  Figure  6 
we  show  model  Z  and  two  additional  models.  CIT11CS3  is  a  smoothed  version  of  CITll. 
Model  Z1  is  a  more  abrupt  version  of  model  Z.  Both  differ  from  model  Z  primarily 
in  the  location  and  abruptness  of  the  discontinuities.  Neither  are  designed 
specifically  to  satisfy  our  present  data  but  merely  to  illustrate  the  difficulty 
of  interpretation  when  discor. tinuities  ar*>  present.  Figure  7  gives  the  computed 
dr /dJjs  tor  these  three  models.  Clearly,  later  arrivals  rust  be  used  to  discuss 


details  of  transition  zones. 
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TABLE  1 


Date 

O.T. 

h:m 

A 

(°) 

Az 

(°) 

h 

(km) 

H 

d  t  /d  A. 

sec/deg 

1963 

Oct  4 

08:34 

22.3 

141 

41 

4.5 

10.9 

7 

05:35 

14.6 

172 

33 

4.7 

13.0 

11 

09:51 

17.0 

161 

33 

4.6 

12.3 

13 

16:45 

14.4 

SE 

50 

4.8 

11.7 

19 

01:04 

24.1 

SSE 

34 

4.3 

8.8 

29 

07:01 

12.3 

NW 

38 

4.7 

13.8 

Nov  14 

09:05 

24.7 

136 

33 

4.9 

9.0 

20 

01:13 

23.5 

330 

39 

4.0 

10.2 

Dec  3 

12:23 

17.0 

160 

33 

4.0 

11.9 

4 

04:18 

15.0 

169 

33 

4.8 

13.8 

12 

00:38 

16.0 

165 

33 

4.5 

13.0 

15 

16:54 

21.7 

142  • 

35 

4.2 

11.2 

22 

01:20 

25.8 

136 

33 

4.4 

9.1 

25 

23:02 

15.5 

170 

33 

3.7 

13.7 

28 

06:57 

26.2 

134 

33 

4.5 

9.4 

29 

04:34 

15.9 

161 

33 

3.8 

11.9 

06:22 

16.3 

161 

33 

4.4 

13.8 

31 

21:41 

20.7 

145 

33 

4.3 

11.1 

1964 

Jan  1 

09:43 

16.7 

162 

33 

4.4 

13.0 

19:45 

12.2 

156 

33 

4.0 

13.8 

6 

23:37 

15.2 

169 

33 

4.1 

13.8 

11 

00:40 

21.3 

144 

33 

4.5 

10.7 

17 

22:18 

20.9 

144 

33 

4.2 

11.1 

18 

14:40 

24.5 

137 

33 

4.6 

9.4 

Date 

O.T. 

4 

Az 

h 

M 

d  f  /d  a 

h;  tn 

(°) 

(°) 

(km) 

sec/deg 

1964 

Jan  20 

20:30 

20.8 

143 

33 

4.3 

11.1 

23 

01:42 

16.1 

167 

33 

3.8 

13.4 

r.8 

04:56 

14.5 

313 

17 

4.5 

14.2 

31 

17:07 

19.6 

334 

14 

4.2 

11.1 

Feb  3 

02:00 

26.0 

137 

34 

4.3 

9.1 

17:34 

20.3 

144 

33 

4.1 

10.8 

8 

09:59 

26.8 

133 

33 

4.1 

8.8 

9 

10:49 

28.9 

134 

33 

3.9 

8.3 

13 

02:22 

26.6 

160 

33 

4.6 

8.9 

02:30 

26.8 

161 

33 

4.7 

8.4 

02:39 

27.4 

162 

33 

4.2 

8.5 

14 

12:07 

14.7 

313 

33 

4.5 

13.5 

18 

09:07 

26.1 

161 

33 

4.0 

9.3 

20 

12:03 

25.7 

136 

49 

4.4 

9.2 

21 

03:15 

25.4 

134 

33 

3.8 

8.8 

26 

20:32 

12.1 

303 

27 

4.6 

15.1 

27 

11:35 

16.6 

155 

33 

4.5 

12.9 

Mar  3 

07:04 

16.3 

311 

33 

4.3 

12.4 

20:02 

12.5 

303 

33 

4.8 

14.5 

5 

11:42 

12.2 

299 

33 

4.3 

14.3 

8 

19:12 

14.2 

158 

33 

4.1 

14.4 

19:44 

14.0 

156 

33 

4.0 

14.3 

10 

09:04 

15.2 

170 

33 

4.3 

13.8 

12 

18:35 

14.3 

171 

33 

3.9 

13.4 

13 

11:54 

28.5 

133 

128 

4.9 

8.8 

Date 

O.T. 

A 

Az 

h 

M 

di  /dA 

h  m 

<°> 

(°) 

(km) 

sec/deg 

1964 

Mar  13 

21:08 

27.0 

132 

33 

4.8 

8.2 

18 

11:39 

26,1 

134 

31 

4.0 

9.4 

23 

06:12 

18.6 

161 

33 

4.0 

11.6 

24 

09:37 

21.4 

327 

22 

4.2 

10.6 

28 

10:08 

11.5 

38 

16 

5.1 

14.0 

April  3 

05:56 

24.4 

136 

33 

4.2 

9.4 

8 

03:38 

23.4 

317 

33 

4.3 

9.8 

June  3 

11:41 

15.0 

169 

33 

4.1 

13.7 

13:50 

13.3 

301 

33 

5.4 

14.7 

4 

04:28 

19.2 

148 

22 

4.7 

10.9 

20 

17:12 

16.5 

160 

28 

5.5 

13.0 

19:35 

16.3 

160 

33 

4.5 

13.0 

24 

15:13 

19.5 

147 

33 

4.4 

11.6 

July  6 

07:22 

18.6 

146 

100 

6.3 

11.3 

14 

12:47 

13.6 

308 

33 

5.4 

14.6 

20 

18:49 

14.6 

172 

33 

5.1 

13.0 

21 

01:07 

14.8 

172 

33 

5.3 

13.4 

03:41 

15.5 

173 

33 

3.7 

13.2 

Aug  13 

06:35 

14.1 

309 

33 

4.9 

14.0 

08:50 

13.7 

310 

33 

4.9 

14.0 

Oct  22 

Salmon 

18.5 

94 

11.9 

TABLE  2 


(kra/sec) 


depth 

Jeffreys 

X 

Gutenberg 

Y 

Z 

(km) 

0 

6.1 

6.1 

6.1 

6.1 

6.1 

30 

6.1 

6.1 

6.1 

6.1 

6.1 

30 

7.75 

8.08 

60 

7.84 

7.87 

100 

7.95 

7,83 

160 

8.13 

8.00 

200 

8.26 

8.12 

225 

8.34 

8.18 

240 

8.38 

8.24 

8.25 

260 

8.44 

8.33 

8.32 

280 

8.52 

8.41 

8.42 

300 

8.58 

8.60 

8.51 

8.50 

8.53 

320 

8.68 

8.59 

8.64 

340 

8.77 

8.68 

8.83 

360 

8.86 

8.78 

8.94 

380 

8.96 

8.88 

9.04 

400 

8.90 

9,05 

9.00 

8.98 

9.13 

420 

9.15 

9.08 

9.21 

440 

9.25 

9.18 

9,29 

4t)0 

9.35 

9.29 

9.38 

480 

9.45 

9.40 

9.46 

500 

9.55 

9.49 

9.50 

9.55 

520 

9.66 

9.61 

9.65 

(km/ sec) 


depth 

Jeffreys 

X 

Gutenberg 

Y 

Z 

Vkm) 

540 

9.91 

9.77 

9.73 

9.74 

560 

9.87 

9.84 

9.84 

580 

9.98 

9.95 

9.94 

600 

10.25 

10.09 

9.99 

10.06 

10.05 

620 

10.20 

10.18 

10.16 

640 

10.31 

10.29 

10.28 

660 

10.42 

10.41 

10.50 

680 

10.54 

10.52 

10.62 

700 

10.65 

10.63 

10.74 

720 

10.76 

10.75 

10.85 

740 

10.88 

10.86 

10.95 

760 

10.99 

10.97 

11.05 

780 

11.10 

11.07 

11.14 

800 

11.00 

11.21 

11.00 

11.18 

11.24 
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FIGURE  CAPTIONS 
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Figure  i. 
Figure  3 

Figure  4 

Figure  5 

Figure  6 
Figure  7 


Location  cf  epicenters  and  array.  Insert  shows  the  geometry 
and  orientation  of  the  array. 

Representative  seismograms. 

Apparent  reciprocal  velocity  versus  distance.  Velocity 
model  Z  is  based  on  the  segmented  curve.  The  dashed 
curve  is  for  the  Jeffreys -Bullen  model. 

Standard  velocity-depth  curves  and  three  structures 
determined  in  this  paper. 

Travel-time  residuals  from  the  J.  B.  tables  for  the  three 
new  velocity  structures. 

Model  Z  and  two  other  hypothetical  mantle  models. 

Theoretical  apparent  reciprocal  velocity  curves  for  the 
models  shown  in  Figure  6. 
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